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ABSTRACT 

In this paper considered the model of two competing populations with double nonlinear diffusion and three types 
of functional dependencies. The first dependence corresponds to Malthusian type of demographic processes, 
the second - Ferkhulst (logistics population), and the third the population of the "Allee" type. The common element of this 
kind of description is the presence of a linear source. And in the descriptions of populations of Ferkhulst and Allee type 
also presents nonlinear sinks. Proposed appropriate initial approximation for quickly convergence iterative process. 
Numerical experiments are carried out with visualization for different values of the system parameters. Modeling of the 
processes of growth of dissipative structures in reaction-diffusion (RD) systems contributes to the development of the 
theoretical ideas about the colonial organization of populations. 

KEYWORDS: Biological Population, Nonlinear System of Differential Equations, Initial Approximation, Numerical, 
Iterative Process, Self-Similar Solutions 

1. INTRODUCTION 

Population model of Kolmogorov Fisher type with nonlinear cross-diffusion considered in [1-12]. 

Let us explain what we mean by the term "cross-diffusion" (or cross-diffusion). Consider the following system of 
two equations in one-dimensional case: 



du x 3 H t ,9 



f a,. > 



= f(u 1 ,u 2 ) + D l —— + h l 



dt 12 1 dx 2 1 dx 



Ql (Ui,u 2 ) — 
ox 



v J (1) 

du 2 d 2 u 2 , 3 ( _ , , du, ^ 



= g(u 1 ,u 2 ) + D 2 — -ir + h 2 — QAu^Uj)— 1 - 



Bt dx dx 



V UA J 



dx 



Cross-diffusion means that spatial move one object, described one of the variables is due to the diffusion of 
another object, described by another variable. At the population level simplest example is a parasite (the first object, 
located within the "host" (the second object) moves through the diffusion of the owner. The term "self-diffusion" 
(diffuse, direct diffusion, ordinary diffusion moves individuals at the expense of the diffusion flow from areas of high 
concentration, particularly in the area of low concentration. The term "cross-diffusion" means moving/thread of one 
species/ substances due to the presence of the gradient other individuals/ substances. The value of a cross-diffusion 
coefficient can be positive, negative or equal to zero. The positive coefficient of cross-diffuse indicates that the movement 
of individuals takes place in the direction of low concentrations of other species occurs in the direction of the high 
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concentration of other types of individuals/ substances. In the nature of the system with cross-diffusion quite common and 
play a significant role especially in biophysical and biomedical systems. 

Equation (1) is a generalization of the simple diffusion model for the logistic model of population growth [13-16] 
of Malthus type (f l (u l ,U 2 } = U l ,f l {u l ,U 2 ) = U 2 ,f 2 (u l ,U 2 ) = U l , f 2 (u 1 ,U 2 ) = U 2 ), Ferkhulst type 
(/[(«,, U 2 ) = U 1 (l-U 2 ),f 1 (u 1 ,U 2 ) = U 2 (l-U 1 ),f 2 (u l ,u 2 )=u i (l-u 2 ),f 2 (u l ,U 2 ) = U 2 (l-U 1 )), and Allee type 
(/i(«i,« 2 ) = u l (l-u/'),f l (u l ,u 2 ) = w 2 (l-w/ 2 ), f 2 (u 1 ,u 2 ) = u l (l-u 2 A ) , f 2 (u x ,u 2 ) = u 2 (\-uf 2 ), 
/?[ > 1, f3 2 >\) for the case of double nonlinear diffusion. In the case when f3 v > 1, J3 2 > 1 , it can be regarded as the 
equation of nonlinear filtration, thermal conductivity at simultaneous influence of the source and the absorption capacity of 

p p, du, du 2 

which is equal to respectively U l ,—U 2 , U 2 ,—U l , exposed by convective transfer rate — — and — — . 

dx dx 

Consider the three-dimensional analogue of the Volterra-Lotka competition with nonlinear power dependence of 
the diffusion coefficient on population density. In the simplest case of Volterra competitive interactions between 
populations can construct numerically, and in some cases analytically inhomogeneous in space solutions [19]. 

2. LOCALIZATION OF THE WAVE SOLUTIONS OF REACTION-DIFFUSION SYSTEMS WITH 
DOUBLE NONLINEARITY 

Consider in the domain Q={(t,x): 0< t < °o, xeR 2 ) parabolic system of two quasilinear equations of reaction- 
diffusion tasks of biological populations of Kolmogorov-Fisher type 



f 



dt dx 



du. 



r* Wi -1 

D x u 2 1 



du x 


p~ 2 3 ^ 
ou x 


dx 


dx 

J 



+ 



dt dx 



2 U \ 



m ? —\ 



du- 



dx 



P-2 \ 

ou 2 
dx 



dx 

+ l(t)^ + k 2 (t)u 2 (l-uf 2 ), 
dx 



(2) 



M l|r=0 M io(-"-)' W 2 1 /=0 M 20 ( X ) ' 

which describes the process of biological populations in nonlinear two-component environment, and its diffusion 



coefficient is equal to D X U 2 



du. 



dx 



p-2 



2 U \ 



m ? -\ 



du* 



dx 



p-2 



_and convective transport with speed l(t) , 



where m l ,m 2 , p, /3 X ,J3 2 - positive real numbers, U x =U l (t,x)>0, U 2 =U 2 (t,x)>0 - required solutions. 

The Cauchy problem and boundary problems for the system (1) in the univariate and multivariate cases 
investigated by many authors [15-21]. 

Purpose of this work is to find estimates for solutions and emerging with a free boundary that gives the chance to 
choose the appropriate initial approximation [15] for each value of the numeric parameters. 

It is known that nonlinear equations have wave solutions in the form of diffusion waves. Under the wave meant 
self-similar solution of equation (2) in the form 
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u(t,x) = /(£), 4 = ct + x 

where the constant C is the wave speed. 

Let's build self-similar system of equations (2) - more simple for research systems of equations. Construct a 
self-similar system of equations by nonlinear splitting method [15]. 

Replacing in (2) 



u 1 (t,x) = e 0 V 1 (T(t),T}),T} = X-j l(£)dC, 



-\k 2 (Od( t 

u 2 (t,x) = e 0 v 2 (t(t),Tj) , J] = x-\ l(C)dg , 

0 

lead (2) to the form: 

m,-l ^V, I ,. /^„[(2-/7)*,+(A-m,+l)* 2 >, , B, 



dv 1 


_ a 




drj 


dv 2 




dr 


d/7 



dv x 


P-2 ^ 


d/7 





D«7 
2 V 1 



dv 0 



P-2 -, A 
dv. 



-fc 2 (Oe [(A - m2+1Ul+(2 -" )ij] 'vf 2 v 2 , 



(3) 



ViL)=v 10 (7), v 2 Lo=v 20 07) 



If 



k x (p - (m x + 1)) = k 2 (p- (m 2 + 1)) . 



then 



by 



r(t) 



[(m, -1)* 2 +(p-2)*, Jf [(m 2 -1)*, +(p-2)£ 2 ]r 



(m, - \)k 2 + (p — T)k x (m 2 - \)k x +(p — 2)k : 



, we get the following system of equations: 



dV[ _ d 
3r drj 



dv 2 _ d 



m, —1 



3vj 


P-2 "N 

dV[ 


8/7 


57 y 



3t drj 



Dm 
2 V 1 



dv,, 



d^ 



p-2 -, A 
dv. 



d^ 



-a, {t)T h v lV ?, 
-a 2 {t)r b2 v^-v 2 , 



where a, =^-2)*, + («, -1)* 2 )\ = (2 "^' +(A + ^ 

(p-2)k x +(m l —l)k 2 

a 2 = kl ( imi -» k ,Hp-2)k 2 r,b 2 = a-*+K+Q-P*> . 

(m 2 - i)k x +(p — 2)k 2 

If b { = 0 , and a ; (0 = const , Z = 1,2 , then the system has the form: 



selecting 



(4) 
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3vj _ d 



dt Btj 
dt drj 



Dm. 



3vj 


p-2 "\ 







Dm 

2*1 



3v 2 


p_2 ^ ^ 


3/; 





-fl 2 vf 2 v 2 . 



With the purpose of obtaining self-similar system for the system of equations (4) firstly we find solution of the 
system of ordinary differential equations 



dv 



dt 
dv 



1 = -a l v l vf 1 , 



= -a 2 vpv 2 , 

dt 
in the form 

V 1 (t) = c l (t + T Q y r ',v 2 (t) = c 2 (t + T Q y^,T 0 >0, 
where 



1 1 1 

1 A A 



And then the solution of system (3)-(4) is sought in the form 



V l (t,Tj) = Vi(t)w l (T,Tj), 
V 2 (t,J]) = V 2 (t)w 2 (T,Jj), 

and X = l(t) is selected as 



r 1 (r) = }v<"- 2 )(Ov 2 <m - 1) (f)* = 



(5) 



1 



:(T + T) 



l-[r,(p-2)+)' 2 (m l -l)] 



l-[y I (/»-2) + y 2 (j»i 1 -l)] 
ln(T + t), if l-i ri (p - 2) + y 2 (jB, - 1)] = 0, 

(T + t), if p = 2 and m,=l, 



! /l-[r 1 (/7-2) + ^ 2 (/n 1 -l)]^0, 



if ^(p - 2) + r 2 K - 1) = r 2 o> - 2) + ^(m 2 - 1) . 

Then for W.(T,x), f = 1,2 we obtain a system of equations 
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dt drj 



dw 2 _ d 



Dm,— I 
1^2 



dw. 



where 



Dm 
2*1 



9wj^ 
drj 



+ ^ 1 (w 1 w 2 A -w,) 



87 



(6) 



drj 



+ y/ 2 (w 2 wf 2 -w 2 ) 



¥i=< 



1 



(l-[y l (p-2) + y 2 (m l -l)])t 

[n(p-2)+r2(»n-l)]\ 



ri«r ( 



^ 2 



1 



(l-[y 2 (p-2)+y l (m 2 -l)])T 

v fr(Hy 2 (p-2)+ri(M2-i)D 
/ 2'n > 



1/ i-[r 1 (p-2)+r 2 (m 1 -i)>o, 

Cri-[r 1 (/»-2) + / 2 (m 1 -l) = 0, 

(/•l-[7 2 (/7-2) + 7 1 (m 2 -l)]>0, 

r i-[r 2 (p-2)+r 1 (w 2 -i)]=o. 



(7) 



Representation of system (2) in the form (5) suggests that, when T — > 00 , Iff • — > 0 and 



3wj 


_ 5 




drj 




_ d 


3r 


drj 



r\ ffk — 1 



dw. 



D 2 w x 



ntj—1 



dw. 



V 



877 



P-2 > 



377^ 

dw 2 
drj 



p-z. 



(8) 



Therefore, the solution of system (1) with condition (5) tends to the solution of the system (8). 
If l — [y 1 (p — 2) + y 2 (m l - 1) = 0 , wave solution of system (6) has the form 

w i (r(t),ri) = y i (Z), g = cT±1], / = 1,2, 

where c - velocity of the wave, and the fact that the equation for w i (T, x) without the younger members always 
has a self-similar solution if 1 — [%( p — 2) + y 2 {m x — 1) ^ 0 we get the system 



d_ 
d% 

d_ 



(v 



m, — 1 



dy x 



m-» — 1 



dy 2 



p-2 



d% 



p-2 



J) + ^ + r ,(y,-y, y «) = 0, 



After integration (8) we get the system of nonlinear differential equations of the first order 
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m,— 1 

y 2 



m 9 — 1 



d£ 

dy 2 



p-2 



d£ 



p-2 



dy 2 



dl; 



+ cy t =0, 



cy 2 = 0. 



(9) 



System (9) has an approximate solution of the form 

y x = A{a-^, y 2 = B(a-^Y\ 
where 



n = 



(p-l)(p-(m 1+ l)) 
(p-2) 2 -(m 1 -l)(m 2 -l) 



(p-D(p-(m 2 +l)) 
72 (p-2) 2 -(m 1 -l)(m 2 -l)' 

And the coefficients A and B are determined from the solution of a system of nonlinear algebraic equations 

(r l y~ l A p - l B m> - 1 = c , 
(r 2 y~ l A m2 - l B p - 1 = c . 

Then taking into account expressions 
u l (t,x) = e° v^tit),?]), 

t 

u 2 (t,x) = e 0 v 2 (t(t),J?) 
we have 



u l (t,x) = Ae" (cz(t)-^l 1 , 



u 2 (t,x) = Be° (ct(t)-^Y + 2 , c>0. 

Due to the fact that 

t 

[bT(t)-^l(7])d?]-x] = 0, 



if 
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x>[bT(t)-$l(Tj)d?]-x]<0, \/t>0, 



then 



u 1 (t,x) = 0, u 2 (t,x) = 0,x>[bT(t)-^l(j])dJ]-x]<0, Vf >0. 

0 

Therefore, the condition of localization of solutions of (2) are the conditions 

e 

jl(y)dy < 0 , Tit) < oo for \/t > 0 . 



(10) 



Condition (10) is the condition for the emergence of a new effect - the localization of the wave solutions (2). 
If the condition (10) unfulfilled, then there is the phenomenon of the finite speed of propagation of a perturbation, i.e. 

{ 

t (m l +p-3)jk 1 (y)dy 

W ; (t, x) = 0 at pc > b(t) , T(t) = j e — d£ , and front goes arbitrarily far away, with increasing time since 

o 

Z(t) — > °° at t — > oo . 

Below we consider the properties of solutions of systems of parabolic equations with cross diffusion. 

3. POPULATION MODELS OF KOLMOGOROV-FISHER TYPE WITH NONLINEAR 
CROSS-DIFFUSION 

Let's consider in the domain Q={(t,x): 0< t < oo, xeR 2 } parabolic system of two quasilinear equations with 
nonlinear cross - diffusion 



dt dx 



du 2 d 
dt dx 



D x u 2 ' 



dx 



P-2 -j \ 



dx 



dx dx 



+ /(0 — + ifc,(0"i(l-"f ) 
dx 

+ l{t) + k 2 {t)u 2 (l - uf 2 ) 
dx 



(11) 



1 /=o~ 



u l0 (x) , u 2 \ t= Q— u 20 (x) . 



which describes the process of biological populations in nonlinear two-component environment, which diffusion 



coefficient is equal D { U 2 



dx 



P-2 



du, 



dx 



P-2 



and convective transfer with speeds l(t) , 



where m l ,m 2 , p, J3 X ,/3 2 - positive real numbers, U l =U l (t,x)>0, U 2 =U 2 (t,x)>0 - search solutions. 
Let's build self-similar system of equations (1 1) - more simple for research of systems of equations. 
Construct a self-similar system of equations by nonlinear splitting method [15]. 
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Replacement in (1 1) 

-J*i(0<*f * 
u x (t, x) = e 0 v, (t(t), Tj),J] = x-j l(C)dC , 

o 



-\k 2 (()di t 

u 2 (t,x) = e 0 v 2 (T(t),Jj),J] = x-j l(C)dC, 



lead (1 1) to the form: 



3V[ _ 3 
dz df] 

3v 2 _ 3 



D,v 2 ' 



3v, 



Dm- 
2V 



3/7 

3v, 



3/7 



3v 2 
3/7 

3v 



J 

1_ 

3/7 y 



■fc 1 (f)e [(A " m '" p+2U2 ^' l 'v 1 v 2 A , 



-k 2 (t)e l(A - m i- p+2)k '- k > ] 'vf 2 v 2 



(12) 



v iL=o =v io( ; 7). v 2 L 0 =v 20 (//-). 



By choosing f(t) = 



[(m 1 +p-2)* 2 +t,]» 



[(m 2 +p-2)t!+i 2 ]J 



(mj + /? - 2)& 2 + fc, (m 2 + /? - 2)^ + k. 



, we obtain the following system of equations: 



3vj _ 3 
3r 3/7 

3v 2 _ = _3_ 
3r 3/7 



r-i /f 

£>,v 2 



3v- 



3/7 

3v, 



3/7 



3v 2 ^ 
3/7 

3vi_ 
3// 



-a, (f)*V* 



(13) 



where a x = fc, ((m l + p — 2)k 2 — k x f 1 , h x — 



(A ~ m i ~ P + 2 )^2 ~ &i 

(m l + p - 2)k 2 - k x 



(m 2 + p — 2)k x —k 2 
If b i = 0 , and a t (t) = const , i = 1,2 , then the system has the form: 



dv x _ d 
dv 2 _ d 



j-v m, -1 



3v n 



Dm 
2^1 



3/7 

3v, 



3/7 



dv 2 ^ 
drj 

-2 A 

drj 



Below we describe one way of obtaining self-similar system for the system (13). It consists in the following. 
We first find the solution of ordinary differential equations 
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dv 



dx 
dv 2 
. dr 

in the form 

V, (?) = Cl (T + T 0 Y* , V 2 (T) = c 2 (t + T o y y \T 0 >0, 
where 



1 1 1 

Ci =i. n=-^-, c 2 =l, y 2 =— 

Pi A 



And then the solution of (12) - (13) searched in the form 



V 1 (t,Tj) = Vl(t)w 1 (T,Tj), 
V 2 (t,7]) = V 2 (t)w 2 (T,7]), 

and X = is selected as 



(14) 



T l (T) = \vl p - 2 \t)v^- l) (t)dt-- 



1 



i-[y 1 (n\ + p-2) + y 2 ] 

ln(T + T), 

(T + T), 



(T + T) 



i-in(»>i+p- 2 )+n] 



»\ ifl-[y 1 {m 1 + p-2) + r 2 ]*0, 

if l-[y 1 (m l + p-2) + y 2 ] = 0, 
if m l + p = 2, 



if Y x (m l + p-\) = Y 2 (m 2 + p-\). 



Then for W t (T, X), 1=1,2 we obtain a system of equations 



3Wj _ 3 

3r 87 



3w 2 _ 3 



n "',-1 9w 2 



3/7 



P dw 2 ^ 



dr drj 



D 2 w[ 



m 2 -l 



' dw. 



J 

- \ 

3 77 J 



+ ^ 2 (w 2 wf 2 -w 2 ) 



(15) 



where 



¥1 = 



if l-[/ l (m 1 + p-2) + r 2 >0, 



l-[r 1 (m 1 + p-2) + y 2 ]T 
rA -(r,(» tl+ ,-2) + r 2 ) j ^ i_ [yi ( ini + p _2) + y2 = 0. 



(16) 
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¥2 



l-[y 2 (m 2 + p-2)+y 1 ]t 



(Mmi+p-v+h) 



if l-[y 2 (m 2 + p-2) + y >0, 
if l-[y 2 (m 2 + p-2) + y=0. 



Presentation system (11) as (14) suggests that when f — > °° , \ff t — > 0 and 



dw l _ d 
3f djj 



dw 2 _ d 



Dy 2 



3t 8/7 



D 2 w™ 2 ~ 1 



drj 



p- 2 \ 
dw 2 

drj 

P-2 -s ^ 

drj 



(17) 



Therefore, the solution of system (11) with condition (14) tends to the solution of the system (17). 
Wave solution of the system (17) has the form 

w i (T(t),rj) = y i (£), g = CT±T}, 1=1,2, 

where c - speed of the wave, and the fact that the equation for w t (T, x) without the younger members always has 
a self-similar solution if 1 — [y l (m l + p — 2) + y 2 ^ 0, we obtain the system 



— w 1 



dy 2 



d 



m-! — 1 



dy l 



p-2 



d{ 



p-2 



After integration (17) we get the system of nonlinear differential equations of the first order 



m, — 1 



y 2 



dy 7 



m-, — 1 



yi 



d% 

dy 1 



p-2 



d£ 



p-2 



dy 2 
d£ 

dy l 
d{ 



+ cy x = 0, 



+ cy 2 =0. 



(18) 



P~ l A- l f} m i + P- 2 



The system (18) has an approximate solution of the form 
y^Aia-ZT, y 2 =B(a-£y\ 

where A and B find from the system of algebraic equations (y 2 ) p 1 A l B 

{yy- l A m ^ p - 2 B- x =c, 

7l l-(2-(m 1 + p))(2-(m 2 + p))' 
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(p-l)(l-(m 2 + p)) 
7l \-{2-{m l + p))(2-{m 2 + p)) 

Then taking into account expressions 

u l (t,x) = e° ^(7(0,77)' 

u 2 (t,x) = e 0 v 2 (t(t),J]) 
we have 

u l (t,x) = Ae 0 (cr(f)-£>? , 
u 2 (t,x) = Be 0 (cT(t)-^Y + 2 , c>0. 

The study of qualitative properties of the system (2) allowed to perform numerical experiment based on the values 
included in the system of numerical parameters. For this purpose, as the initial approximation was used asymptotic 
solutions. For the numerical solving of the task for the linearization of system (2) has been used linearization methods of 
Newton and Picard. To build self-similar system of equations of biological populations used the method of nonlinear 
splitting [16,19]. 

4. NUMERICAL EXPERIMENT 

For the numerical solving of the task (2) construct a uniform grid 

0) h — {x; — ih, h>0, i- 0,1,..., n, hn - /}, 
and the temporary net 

°>h , = fr; = A > K > °> J = QX-, n, vn = T\. 

Replace the task (2) by implicit difference scheme and receive differential task with the error 

As you know, the main problem for the numerical solution of nonlinear task is the appropriate choice of the initial 
approximation and the method of linearization of the system (2). 

Consider the function: 

v w (t,x) = v l (t)-(a-^)l 1 , 
v 20 (t,x)=v 2 (t)-{a-g)l 2 , 

where V 1 (t) = e k 'V l {t) h V 2 {t) — e kt V 2 (t) defined above functions, 
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Note (a) + means that (a) + =max(0,6[). These functions are finite speed of propagation of perturbations 
[16,19]. Therefore, for the numerical solving of the task (l)-(2) fS x > <7[ _as an initial approximation offered functions 

v i0 (t,x) ,i = l,2. 

Established in the input language MathCad program allows visually trace the evolution of the process for different 
values of parameters and data. 

Numerical calculations show that in the case of arbitrary values (7 > 0, f5 > 0 qualitative properties of solutions 
are not changed. Below listed results of numerical experiments for different values of parameters (Figure 1 - Figure 3). 



Paramtftr: Value* 


E*mlts ofNumtrical Esptrinent 


M = 2.1 , jHj = 2.3 r p =3 
t =0.1 

epr=10" J 






iLpfj 


to ii J*> 
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M" 
iii* 
it- 
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* ii » ii 



Figure 1: Dynamics of Malthusian Population 



Param*t*r Values 


Re Hilti D-f ^Tiunericil Eiptrimttit 


ff\ =21 .m* =2.J : p=J 
fi=l =7 
A =L k, = ? 
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fc f- 

L 1 














m = 2.3^=2.1^=3.} 
A = 1 fc = 3 




l. 




a: 


\ ,. % ■ 









Figure 2: Dynamics of Logistic Population 
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Parameter Value : 
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Figure 3: Dynamics of Allee Effect 



5. CONCLUSIONS 

In the future will be explored theoretical properties of systems with cross-diffusion, which will make a significant 
contribution to the study of nonlinear systems of differential equations. 
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